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Abstract 

By studying its subgraphs, it is argued that the lower critical percolation 
threshold of the enhanced binary tree (EBT) is bounded as p c i < 0.355059, 
while the upper threshold is bounded both from above and below by 1/2 
according to renormalization-group arguments. We also review a correlation 
analysis in an earlier work, which claimed a significantly higher estimate of 
p C 2 than 1/2, to show that this analysis in fact gives a consistent result with 
this bound. Our result confirms that the duality relation between critical 
thresholds does not hold exactly for the EBT and its dual, possibly due to 
the lack of transitivity. 
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1. Introduction 

Percolation has been always of physical interest since it was introduced 
as a description for a fluid in random media jlj]. The percolation problem 
has served as one of the basic models in understanding critical phenomena 
from a geometric viewpoint jij. An interesting aspect of the percolation 
phenomenon is that a lattice with a constant negative curvature may have 
two different thresholds, which coalesce into one as the curvature vanishes [3[. 
Figure []Ja) gives the simplest example of a hyperbolic lattice: a regular 
tree with level L, where L is defined as the maximum path length from the 
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origin. The set of points at the same distance I from the origin constitutes 
the Zth layer. In this tree, every vertex has three neighbors except at the 
boundary, and the number of vertices scales as iV ~ 2 L . An immediate 
consequence is that the number of boundary points, which scales as B ~ 2 L_1 , 
always occupies a finite fraction of N even in the thermodynamic limit. The 
existence of two percolation thresholds has been attributed to this peculiar 
structural property. That is, when the occupation probability p reaches the 
lower threshold, p&, a system- wide connection is first achieved, while the 
largest cluster occupies a finite fraction of N only if p exceeds the upper 
threshold, p c2 . Therefore, if we count the boundary points connected to the 
midpoint of the system as we increase p, the connection simply does not exist 
at p < Pd, and even after passing this first threshold the connected boundary 
points remain as a negligible part of the whole boundary until p reaches p c2 . 
In the case of the bond-percolation problem in a simple tree as in Fig. QJa), 
for example, one can easily show that p*^ ee = 1/2 and p^ 2 e = 1 0. 

The enhanced binary tree (EBT) is a nontrivial model derived from a 
tree pjj: it is obtained by adding bonds to the tree between every pair of 
neighboring points within each layer / > 1 [Fig. E3Jb)], so it has p c \ < p^ e 
and p C 2 < p^ 2 e due to the existence of loops. An interesting analogy of 
an EBT would be a biological taxonomy which is mostly a tree structure 
but with genes horizontally transferred as well (see, e.g., M). Although the 
lower threshold could be easily measured as p& = 0.304(1) [4|, there remains 
a controversy in locating p c2 

BBS: 1 

claimed that p c2 = 0.564(1), and 
the correlation analysis was suggested to support this claim [7(, while other 
numerical methods preferred p c2 = 0.48(1) Q. Recently, there appeared 
an analytic calculation [8|, which sug gest s a possible generalization of the 
triangle-triangle transformation 0, [lO, 11]: in order to use this transforma- 
tion in an exact fashion, one should be able to decompose a given structure 
of identical triangular unit cells, where the three vertices of one unit cell are 
denoted as A, B, and C, respectively (see Fig. [2](a)). Then P(A, B, C) is 
defined as the probability that A, B, and C are all connected, P(A, B, C) 
as the probability that B and C are connected but A is disconnected from 
them, and P(A, B, C) as the probability that none of them are connected 
to each other. One may locate the critical threshold by equating these two 
probabilities (lb1 |: 

P(A,B,C) = P(A,B,C). (1) 
From a trivial equality, P(A, B, C) + P(A, B, C) + P(A, B, C)+P(A,B,C) + 
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Figure 1: Schematic plots of hyperbolic structures drawn on the Poincare disk, (a) A tree 
structure up to level L = 10 and (b) the enhanced tree obtained from (a). The midpoint 
is located at the origin in each plot. 




Figure 2: (a) A triangular cell having vertices A, B, and C. (b) An array of such triangular 
cells, where the Ti indicate the points in the upper layer on the right-hand side of S. 



P(A, B,C) = 1, one can express Eq. ([I]) as 



P(A,B) P(A,C) 

+ 7 i n\ = * = 1, (2) 



l-P(A,B,C) 1-P(A,B,C) 

where P(a, 0) means the probability of connection between vertices a and 
0. If one picks up an arbitrary starting point S in a very large layer as in 
Fig. Hfb), the probability of connecting to the upper layer on its right-hand 
side can be written as Yl^o Prob(Tj is the first right-hand connection) = 
Y^ =0 P{A,B)P i {A,B,C), which is identical to the first term in Eq. Q2). 
Likewise, the second term in Eq. ([2]) means the probability of a left-hand 
connection from S. In this sense, Eq. ([I]) can be regarded as describing a 
certain connective property between two adjacent layers at criticality (see 
also [12], where one finds a similar idea). Applying this idea to the EBT 
without requiring the self-duality of the triangular type, one obtains ^f(p) = 
p(l+p)/[l—p(l—p)] [8j. Then, the equation ^f(p) = 1 is satisfied at p = 1/2, 
which in [8[ was interpreted as an exact value of p& for the EBT. Since the 
simple binary tree has \I/ tree (p* 1 j ee ) = 1/2, assuming this again holds for the 
EBT, one obtains the value of p cl as (vT3 - 3)/2 w 0.302776 Q. 
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Figure 3: Subgraphs of the EBT for bounding p c \. (a) The tree made of triangles has 
its lower threshold at p\ « 0.403032, which bounds p c % from above, (b) By adding more 
bonds at every second layer, a larger unit cell is constructed, and we get a sharper bound, 
p* 2 ~ 0.373897. (c) The largest unit cell considered in this work, giving the upper bound 
Pel < P3 ~ 0.355059. The dotted lines mean that other unit cells are attached there. 

In this work, we use alternative approaches and argue that p c \ < 0.355059 
by examining solvable subgraphs and p C 2 < 1/2 by means of the renormalization- 



group (RG) method on hierarchical structures [13] . We also demonstrate that 
the correlation analysis indeed gives a consistent result with this argument 
for p c2 . The paper is organized as follows. We explain the results for p cl and 
p C 2 in Sec. [2] and Sec. [3J respectively. Then Sec. H] additionally discusses the 
correlation behavior, and Sec. [5] summarizes this work. 

2. Lower threshold 

A critical percolation threshold becomes higher when part of the links are 
removed from the original graph. We could therefore argue that p c \ < 1/2 
above by considering the simple binary tree as a subgraph of the EBT. We 
refine this bound by taking larger but still solvable subgraphs. 

One may first consider a tree of triangles, which is also a subgraph of the 
EBT (see Fig. 13(a)). We are going to ask whether there can be any path 
penetrating from the top to the bottom. A relevant quantity would then be 
how many vertices can be reached on the (/ + l)th layer from a single vertex 
on the /th layer. If we focus on a unit cell, which is simply a triangle here, we 
find 2 3 possible configurations, since it has three bonds. By checking all of 
these configurations, one can easily get the average number of descendants, 
i.e., the expected number of bottom points connected to the top vertex within 
the cell. It is easily given as 



mip) = 2p(l -p) 2 + 6p 2 (l -p) + 2p 3 = 2p[l +p(l -p)) . 
Solving m(pj) = 1, we obtain p\ = § (l + 2cos# - 2V3sin0) « 0.403032, 
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where 6 = ~ arctan |^^=-J . This provides an improved upper bound for p cl . 

Note also that ni(p) is a monotonically increasing function of p from zero 
to 2. This implies that the upper threshold of this subgraph is located at 
p = 1, since otherwise the fraction of connected vertices will decrease as the 
number of passing layers increases. 

This sort of construction requires that a unit cell should possess only 
one top vertex and that every bottom-layer point equally becomes a new 
top vertex for a subsequent unit cell. We can therefore extend this idea a 
little further, as follows. For every second layer, we add a bond between 
two daughter triangles under the same mother triangle, as in Fig. EJ^b). This 
creates a new unit cell containing three triangles plus one bond between 
the daughter cells. Such a cell has 10 bonds in total, meaning 2 10 = 1024 
possible configurations. Directly enumerating them shows again how the top 
vertex connects to the bottom layer within this cell. The average number of 
descendants then reads 

n 2 {p) = 2p 2 [1 + p(l - p)] [2 + p(l - p)(3 + 3p + p 2 - I0p 3 + 5p 4 )] , 

which ranges over [0,2 2 ]. Numerically solving ^(pg) — 1 leads to p* 2 ~ 
0.373897. It is still possible to consider a larger unit cell with seven triangles 
(see Fig. [3^c)) and find the average number of connected bottom-layer points 

as 

n 3 (p) = 2p 3 [1 + p{l - p)} [4 + p(l-p)(13 + 26p + 23p 2 -37p 3 -64p 4 
-42p 5 - 116p 6 + 234/ + 469p 8 - 830p 9 + 181 lp 10 - 2898p n 
-4735p 12 + 21801p 13 - 31538p 14 + 24399p 15 - 10894p 16 
+2664p 17 - 278p 18 )] 

from 2 25 configurations. By solving ns{p^) = 1, we find an upper bound for 
Pd as p* 3 ~ 0.355059. On the other hand, it is readily seen that p c \ can never 
be lower than 0.25, since every vertex in the EBT has z = 5 neighbors except 
the zeroth layer, so that the tree approximation yields \f(z— 1) = 1/4. We 
note that both the numerically obtained threshold p c \ = 0.304(1) [4j and the 
analytic prediction p c % = — 3)/2 w 0.302776 [8j lie within these upper 
and lower bounds. Since this subgraph method soon becomes impractical due 
to the huge number of possible configurations as the size of a unit cell grows, 
there is no sharper bound available at present. One may try an extrapolation 
by using the correlation- length scaling exponent v — 1 at p — p c \ [3j], but this 
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Figure 4: (a) A part of an EBT structure. The tree part (solid) is occupied with probability 
p and the rest of the bonds, i.e., the horizontal lines, are occupied with probability q. At 
p = 1/2, a single bond bi will directly connect two groups of boundary points, G\ and 
G2, with probability q, while 62 will connect G3 and G4. (b) Schematic description for 
r n +i, the global-connection probability for the bottom layer at the (n+ l)th RG iteration. 
Hence the probability r n at the nth iteration is related to r„ + i by Eq. ([3]). (c) Dynamics 
of the RG flow for q < 1/2 according to Eq. ([3]), which is here drawn as the solid line, 
(d) An alternative recursion scheme illustrating Eq. ((4|), where the double lines represent 
coarse-grained links, i.e., z n or z n+ \. 

will be a rather crude estimation. It is also clear that these subgraphs do not 
delimit p C 2 of the EBT at all since all of them have their upper thresholds at 
p — 1. 

3. Upper threshold 

The main idea in bounding p C 2 begins with the fact that a simple binary 
tree has a lower threshold p* r 1 ee = 1/2. We assign occupation probabilities 
p and q to the tree part and the rest of the EBT, i.e., the horizontal bonds 
in Fig. IHa), respectively. Let us start from q = 0, where the EBT becomes 
identical to a simple tree. At p = 1/2, the tree part reaches its first critical 
threshold so that each vertex is connected to one boundary point on average. 
Then we increase q so that a finite fraction of boundary points can merge 
into a single cluster. 
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An interesting property of the EBT is that the horizontal bonds would 
have different connection ranges depending on which layer they belong to. 
In Fig. II|a), bond b\ directly connects two groups of boundary points, G\ 
and G2, at average distance 2, which is measured along the bottom layer. 
On the other hand, the average range of connection between groups G3 and 
G4 by bond 62 is twice as long as than that. In other words, if we focus only 
on the bottom layer by setting p = 1/2, the horizontal bonds constitute a 
hierarchical structure so that the range of connections becomes twice as long 
every time when a new higher level is introduced. Such a property allows 



one to formulate an RG equation 13j, which here can be written as 



n+l 



q + (l-qK, (3) 



where r n is the probability of global connection in the bottom layer at the 
nth RG scaling transformation (see Fig. H](b)). Here, we are asking ourselves 
how probable it is for two different points that we have arbitrarily picked 
up from the boundary bottom layer to belong to the same cluster. Since a 
cluster containing the middle part of the system occupies only a negligible 
fraction of the whole boundary at p < p c2 by definition, this means that such 
a chance gets significantly large above p C 2, and thereby Eq. ([3]) is related to 
the upper threshold. It is stated in Eq. Q that the global connection can 
be established either by a new long link with probability q or by existing 
shorter links. In the limit of large n, we may set r n = r n+1 = r ro and then 
Eq. ([3]) is easily solved to yield a nontrivial stable solution, r ro — q/(l — q), 
as shown in Fig. H](c) [13]. Noting that r^, is responsible for connecting a 
number of the boundary points to one another at criticality, we find that 
q = 1/2 should be a transition point provided that p is fixed at p^ ee = 1/2. 
Since p and q happen to have the same value here, even if considering the 
homogeneous case where q is always set equal to p, we can conclude that 
p = 1/2 is high enough to connect a significant fraction of boundary points. 
In short, p — 1/2 should be higher than or equal to p C 2, the upper critical 
percolation threshold of the EBT. 

We note that Eq. ([3]) is an approximate description with a coarse-grained 
variable r n , chosen for ease of explanation. It is possible to make an alter- 
native recursion scheme which instead yields a lower bound. The scheme is 
illustrated in Fig. H](d), where q is assumed to be equal to p. There are two 
outer points on the left-hand side and two others on the right-hand side. In 
merging a set of bonds into a single one with connection probability z n+ i, we 
are interested in linking any of the left outer points to any of the right outer 
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points. Obviously, a part of contribution to z n+ i comes from filling the upper 
bond with probability p. Even if it is not filled with 1—p, there remain a few 
more possible cases. Suppose that there happen to be two filled bonds which 
are the right one of the left triangle and the left one of the right triangle, for 
example. It is then enough to have only one z n in between. Likewise, we can 
consider all the other cases and arrive at the following recursion relation: 

z n+1 =p + (l- p )[(l- pfzl + 2p(l - p)z 2 n + p 2 z n ] . (4) 

Solving this by setting z n = z n+ i = Zqo as above, one finds that 

a/1 + Ap - 6p 2 + p A 1+p 
Z °° ~ 2(1 -pf 2(1 - P y 

Again, we see that — 1 at p > 1/2. However, since the renormaliza- 
tion includes connections between the lower outer points even if there is no 
interlayer connection, the result yields a lower bound of the threshold. 

4. Correlation Analysis 

A recent numerical estimation in [3] suggested p C 2 > 0.56, which certainly 
exceeds 1/2 given above as an upper bound. We present a brief discussion 
on this discrepancy. Following [7], let us consider the probability C(l) that 
points at level / in the EBT belong to the same cluster as the midpoint does. 
Such a correlation will be a monotonically decreasing function of /, and it 
converges to a constant if p > p c2 . If we are to take only robust behavior 
insensitive to any particular L, a possible way would be to compare two 
different system sizes, say L = 23 and 24, and take only data points up to 
I = (ma*) where they cease to overlap within error bars. Furthermore, we 
need to exclude I < since the EBT has a lower number of connections 
close to the midpoint, so C(l) may decrease anomalously at small I. For 
the data presented in Fig. [5J we set Z min = 2 for p < 0.55 and Z min = 1 for 
p > 0.55 because of the more rapid exponential decay, which possibly implies 
that £ < 2. We assume that C(l) within the range of l m i n < I < Z max will 
show a simple decaying behavior with a certain characteristic length scale £ 
as 

C(l) ~ aexp(-//0 + c, (5) 

where a and c are independent of I. The parameter c can be also said to 
be the asymptotic value of correlation at large I. Since this term in the 
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Figure 5: (a) Correlation function measured for L = 24 over 10 4 samples. The occupation 
probabilities are p — 0.49, 0.51, 0.53, 0.55, and 0.57, from bottom to top. The solid curves 
are according to Eq. ([5]) with least-square-fitted parameter values, while the dotted curves 
are found when we fit the data without c. (b) Correlation length £ obtained from the 
measurements. One observes a peak around p = 0.51 from the fitting by Eq. J5} (solid). 
If one fits the same data after setting c = 0, £ can be overestimated in the supercritical 
phase (dotted), (c) Asymptotic correlations obtained by Eq. ([5]). (d) Reduced \ 2 values 
for c and c = 0, respectively. 
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fitting procedure was not included in [7J, £ is presumably overestimated in the 
supercritical phase, since the nonvanishing part of C (I) would be interpreted 
as a very slow decrease (see Fig. Ufa)). Our fitting results plotted in Fig. [5(b) 
and Fig. [5]^c) show signatures of the transition around p = 0.51, which is 
significantly lower than p = 0.56 claimed in and is in fact fairly close to 
1/2. In order to compare the goodness of fits, we calculate the x 2 statistic 
defined as 



where Oj is the zth data point observed with variance of, and Ei is the 
corresponding expected value from Eq. (jHJ). The number of degrees of free- 
dom, v, is here given as the number of data points minus the number of 
fitting parameters. The reduced \ 2 statistic, x 2 ' l v t allows one to compare 
the performance of different fitting functions, and a rule of thumb states 
that a good fit is achieved when x 2 / u ~ 0(1). In Fig. [3(d), we find that 
including the asymptotic correlation term, c, indeed describes the behavior of 
C (I) better, since it makes the reduced x 2 statistic maintain its value around 
0(1) throughout the checked range of p, while the pure exponential function 
without c becomes a poor description for the same data at p > 0.54. One 
therefore finds that the threshold value determined numerically is consistent 
with the analytic prediction given above, although the result is subject to a 
larger uncertainty than the case of p c \. 

5. Summary 

In summary, we set upper bounds for p cl and p c2 of the EBT structure; 
that is, Pd < pi = 0.355059 and p c2 < 1/2. In addition, we obtained 
Pc2 > 1/2 as well, which confirms the argument in that p c2 = 1/2. Since 
the upper threshold has been particularly under debate, we tried to settle 
the issue to a large extent by showing that the RG method developed for hi- 
erarchical structures led to the bounds of p c2 . We also demonstrated that the 
correlation analysis yielded a consistent result with this RG argument. Both 
of these analytical and numerical approaches disprove the duality conjecture 
that the EBT and its dual lattice are related in such a way that the upper 
threshold for the EBT, p^. and the lower threshold for the dual lattice, p c i, 
sum up to p c2 + pd = 1 j4J, since the lower threshold for the dual lattice is 
numerically determined as p c \ = 0.436(1) to good precision Q. The possible 
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value of p c2 given in this work clearly shows that p c i +p C 2 < 1, as already re- 
ported for other hyperbolic lattices j^]. This observation confirms that such 
a duality relation requires transitivity 14j, which does not hold for the EBT. 

From a methodological viewpoint, our analysis in this work relies largely 
upon the regularity of the underlying structure. There also exist other regular 
hierarchical structures such as the Apollonian networks 
networks 



15, 16 and flower 
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where the percolation problem has been studied by means of 
recursion. The extension and usefulness of our subgraph analysis for these 
cases remain to be investigated. 
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